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In this paper, we investigate the relations between global and local energy transfers in a turbulent 
von Karman flow. The goal is to understand how and where energy is dissipated in such a flow and 
to reconstruct the energy cycle in an experimental device where local as well as global quantities can 
be measured. In order to do so, we use PIV measurements and we model the Reynolds stress tensor 
to take subgrid scales into account. This procedure involves a free parameter that is calibrated using 
angular momentum balance. We then estimate the local and global mean injected and dissipated 
power for several types of impellers, for various Reynolds numbers and for various flow topologies. 
These PIV-estimates are then compared with direct injected power estimates provided by torque 
measurements at the impellers. The agreement between PIV-estimates and direct measurements 
depends on the flow topology. In symmetric situations, we are able to capture up to 90% of the 
actual global energy dissipation rate. However, our results become increasingly inaccurate as the 
shear layer responsible for most of the dissipation approaches one of the impellers, and cannot be 
resolved by our PIV set-up. Finally, we show that a very good agreement between PIV-estimates 
and direct measurements is obtained using a new method based on the work of Duchon and Robert 
[J. Duchon and R. Robert, Nonlinearity, 13, 249 (2000)] which generalizes the Karman-Howarth 
equation to nonisotropic, nonhomogeneous flows. This method provides parameter-free estimates 
of the energy dissipation rate as long as the smallest resolved scale lies in the inertial range. These 
results are used to evidence a well-defined stationary energy cycle within the flow in which most of 
the energy is injected at the top and bottom impellers, and dissipated within the shear layer. The 
influence of the mean flow geometry and the Reynolds number on this energy cycle is studied for a 
wide range of parameters. 


I. INTRODUCTION 

Understanding how and where energy is dissipated in turbulent flows has been a great challenge for many years, 
and would have important implications in many areas such as fundamental research, aeronautics or industry. In the 
classical three-dimensional turbulence phenomenology, energy is injected at large scales by the forcing mechanism, 
transferred downscale at a constant rate e following a self-similar cascade, and then dissipated into heat at the 
Kolmogorov length scale, where viscous effects become dominant. In ideal stationary, homogeneous and isotropic 
turbulence, the measurement of energy dissipation can therefore be achieved via 3 independent and equivalent means: 
i) by monitoring the injected energy; ii) by monitoring the dissipated heat; iii) by monitoring the cascade energy rate 
via multi-scale single points measurements of the velocity (via e.g. anemometers or array of hot wires). The first 
two measurements are global, the last one local, but given the homogeneity in space and time, they all provide the 
same information. In most realistic situations, however, the turbulence is anisotropic and/or inhomogeneous and/or 
non-stationary. In such cases, there is not any necessary equivalence between global energy injection, global energy 
dissipation and local energy dissipation. The study of these three quantities requires detailed knowledge of the forcing, 
the heat distribution and velocity over the whole domain. 

In that respect, it is interesting to focus on intermediate situations, where the turbulence is generated by a well 
controlled mechanism, in a simplified geometry, as achieved for example in classical laboratory experiments such as 
Taylor-Couette or von Karman set-ups. In such cases, it is easy to monitor the forcing and implement a cooling 
mechanism so as to achieve stationarity, where global energy injection and dissipation equilibrate on average. The 
local energy dissipation can then be computed from stationary energy budgets derived from Navier-Stokes equation, 
using local measurements of velocity obtained for example using the now classical Particle Image Velocimetry (PIV). 
This technique provides measurements of the instantaneous velocity field Ui at several points of a plane (or of a 
volume) at the same time. From this, one may compute the Reynolds stress tensor Sij = ^{diUj djUi), and study 
the local dissipated power = 2vSijSij. A detailed comparison between this local estimate and the global energy 
dissipation rate in a Taylor-Couette flow at various Reynolds number has recently been made by Tokgoz et al. [T] 
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using tomographic PIV measurements. They observe that the local and global estimates coincide within 10% as long 
as the laminar flow and Taylor vortex flow regimes are fully resolved. However, as the Reynolds number increases, 
the dissipative scale decreases and becomes smaller than the flnite resolution of the PIV (set by the camera resolution 
and the velocity reconstruction algorithm). In that case, the estimate of local energy dissipation based on the velocity 
gradients and on the viscosity ty becomes increasingly inaccurate, and underestimate the global energy dissipation. In 
order to remedy this problem, it has been suggested to use techniques borrowed from Large Eddy Simulation (LES) 
[5] . This allows us to model the subgrid scales (SGS) in terms of the large scale velocity field resolved by the PIV and 
allows computation of all terms in the energy budget Hi] including the terms responsible for the scale-to-scale energy 
transfer. By the Richardson-Kolmogorov cascade picture, this allows an estimate of the local energy dissipation as 
long as the scale used in the computation lies in the inertial range. 

In this paper, we test these methods in a turbulent flow generated by two contra-rotating impellers (von Karman flow) 
for Reynolds numbers ranging from 10^ up to more than 10®. At such Reynolds numbers, the dissipative scale ranges 
from a few millimeters to a few tens of microns. With fixed velocity of the impellers and cooling, the resulting flow is 
stationary, highly anisotropic and inhomogeneous [7] , thereby providing a unique laboratory flow to test local energy 
dissipation procedure. This closed flow geometry permits direct estimates of the global energy injection by torque 
monitoring at the two impellers. Using a large scale Helium facility with calorimetric measurements [S]) were able 
to show that in a stationary state, the global injected power and the global dissipated heat coincide within a few 
percent at large Reynolds number, and for a wide range of operating conditions (including differential rotation of the 
impellers). Moreover, the simple cylindrical geometry allows for Stereoscopic PIV measurements of the velocity field 
over a vertical plane spanning the whole experiment, at a resolution of a few millimeters. Since measurements have 
only been made in a meridional plane, we do not have access to orthoradial derivatives. Hence, it is challenging to 
have an accurate estimate of the local energy dissipation rate despite this lack of information. To do this we first use 
a LES method using statistical axisymmetry and angular momentum budget to calibrate the model. Then, we test 
a generalization of the Karman-Howarth formula derived by Duchon and Robert H that provides a parameter-free 
estimate of the local energy dissipation for any nonhomogeneous, nonisotropic flow. Both estimates are then averaged 
over the whole volume for comparison with the global estimate of the energy dissipation based on torque monitoring. 
This is done for a wide range of control parameters, varying both the Reynolds number, the mean flow geometry and 
the flow asymmetry. We then use these measurements to evidence a stationary energy cycle within the flow where 
energy is injected at the top and bottom impellers, and dissipated within the shear layer. 

This paper is organized as follows: in Section H, we summarize the theoretical tools needed for the implementation 
of the LES technique and the Duchon-Robert (DR) formula in our analysis of the energy dissipation rate. In Section 
HI, we review the von Karman geometry and specialize these formulae to the case of cylindrical geometry. The flow 
diagnostics based on PIV measurements are derived and summarized. In Section IV, we apply these results to a set of 
measurements drawn from our data base of von Karman flow. We first tune our LES model using angular momentum 
balance. Then, we compare PIV estimates of the global dissipated power with direct torque measurements using both 
the LES technique and the DR formula. Finally, in section jV] we evidence the energy cycle of the von Karman flow, 
and we discuss its evolution, as well the evolution of our diagnostics, as a function of the flow topology. 


II. THEORETICAL BACKGROUND AND METHODOLOGY 

In this section, we use a Altering approach to derive energy and angular momentum balance at a given scale ^ from 
Navier-Stokes equations, and we give the expression of the different terms that we will use in our analysis of the 
injected and dissipated power. We consider Navier-Stokes equations 


r dtUi + UjdjUi = —diP + vdjdjUi 

\ djUj = 0 , 

where we use Einstein summation convention over repeated indices. 


( 1 ) 

( 2 ) 


A. The filtering approach 


Following the procedure in cni, we define a coarse-grained velocity field at scale £ as 
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ui{x) = J dr Ge{f)u^{x + r), (3) 

where G is a smooth filtering function, nonnegative, spatially localized and such that J dr G{r) = 1. The function Gi 
is rescaled with I as Ge{r) = £~^G{f/£). Coarse-graining the Navier-Stokes equations gives 

J dtu^ + UjdjUi = — diP^ + vdjdjul (4) 

I d,u] = 0. (5) 

In equation Q we introduced = {uiUjY — ufuj which is the stress tensor from the SGS. 

In what follows, in order to simplify the formulae and for readability considerations, we will drop the index £. Unless 
specified otherwise Ui now denotes the component of the coarse-grained velocity field (same thing for P). 


B. The energy balance equation 

We now take the scalar product of equation Q with Ui, and after a few lines of algebra we get the local energy balance 
equation 


dtE + dj{ujE) = —dj{uiP^) + Sijpi — di{uiP) + v{djdjE + didj{uiUj)) — 2vSijS'‘P (6) 

where E = is the large scale kinetic energy per unit mass, and Sij = ^{diUj + djUi) is the large scale strain rate 
tensor. The overline denotes the statistical average of a quantity. Let us now define a vector = Uj {E-\- P)-\- UiP^ — 
v(djE + di{uiUj)). We can then rewrite (j^ as 


dtE + djJi = Sijpi - (7) 

This local equation is valid in any geometry, for any type of flow, for any filtering with the properties given above. 
In section III D, we apply this formula to the specific axisymmetric von Karman geometry to derive the local energy 
production and dissipation rate per unit mass. 


C. The angular momentum balance equation 

In a very similar fashion, it is also possible to derive an angular momentum balance equation at scale £. We take the 
cross product of rj and equation Q, and after a few lines of algebra we get 


dfEi T — ^ijkdk{rj P) Ctj^dmir + vi^djdjLt + ^ijkdmdkir j'^raY ^^ijkdjUk: (8) 

where etjk is the three dimensional, total antisymmetric Levi-Civita symbol so that Li = CijkrjUk is the component 
of the angular momentum per unit mass. Of course, it can be checked that it is also possible to derive ^ from ^ 
using ^ and the identity E = ^{L^ + (r.uY). As before, we can write (j^ as an equation of conservation for each 
component of L 


dtL, + = 0, (9) 

where = ujL^ - CijkVkP + eimkrmP^ - - e^jkdm{rkUm) + f-ijkUk)- 

In section IV A, we use this angular momentum balance as a constraint to calibrate the LES model, that we now 
describe. 
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D. LES method for balance equations 

The computation of the different terms in Q and ([^ requires the knowledge of both the velocity field at scale I and 
the Reynolds stress In most practical situations, e.g. when the flow is turbulent and the velocity is measured 
through a PIV system, only the former is available, since we cannot resolve the dissipative scale. A traditional way 
to overcome this problem is to use a LES-PIV approach mm to model r*-! in terms of the large scale velocity field 
[in eg. Several models exist. In the present paper, we choose the gradient model m where 


= CAldkUidkUj, ( 10 ) 

where C is a constant to be calibrated and is the width of the filtering. This model ensures forward scatter and 
backscatter of energy between resolved scales and SGS. 


E. Duchon-Robert energy balance equation 

An alternative local energy balance equation has been derived by Duchon and Robert |9] using Leray’s weak solution 
formalism and Onsager’s ideas m- The latter amounts to consider a sequence of coarse-grained solutions of the 
Navier-Stokes equations Eq. Q in the limit i ^ 0 and derive the corresponding energy balance, that reads 

dtE + dj {u^{E + P)- vdjE) = -udjUid^u^ - (11) 

where is expressed in terms of velocity increments 6u(f) = u{x -\- r) — u{x) as 


^(u) = lim ^f(u) = lim — 


dr {'VGe){r) ■ Su(f) \Su{r)\^ 


( 12 ) 


where the dependence of Su and ^ in x is implied. As the Reynolds number tends to infinity, the scale i can be 
chosen as small as one wants, and the quantity can be seen as the contribution to dissipation coming from a 

generalized cascade process (possibly linked with the formation of small-scale singularities). Since the result cannot 
depend on the filtering function Gi, Duchon and Robert specialized the expression to a radially symmetric filter to 
get the alternative expression, devoid of any free parameter 


“ f l<5w(ex)P Su{ex) ■ X , (13) 

IOtt e->o e J\x\=i 

where dE denotes the area measure on the sphere. As noticed by Duchon and Robert, this expression coincides 
with the statistical mean rate of inertial energy dissipated per unit mass derived from the anisotropic version of the 
Karman-Howarth equation. They therefore argue that the previous formula provides a local non-random form of the 
Karman-Howarth equation, valid even for anisotropic, inhomogeneous flows. In the sequel, we apply this formula to 
our PIV measurements to test whether it can provide a parameter-free estimate of the global dissipation, as well as 
local instantaneous maps of the local energy dissipation. 


III. APPLICATION TO A VON KARMAN GEOMETRY 

The goal of the present paper is to compare global estimates of dissipated power with torque measurements, check their 
coincidence, and study maps of local energy dissipation rate to get some insight of the detailed processes governing 
energy transfers in a von Karman flow. 


A. Experimental set-up 

The von Karman experiment has been extensively studied over the past years [ZlIISHSl. We give here a brief review 
of the main features of the set-up. 
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FIG. 1: Illustration of a TP87 type impeller, with 8 blades. The convention used to name the two different forcing conditions 
is represented: (+) when the convex face of the blade goes forward, (-) for the other. The angle a characterizes the curvature 
of the blades and is equal to q = 72° in the case of TP87 impellers. TM60 impellers look essentially the same, with 16 blades 
instead of 8. 


Our von Karman flow is generated in a vertical cylinder by two coaxial, contra-rotating impellers providing energy 
and momentum flux at the upper and lower end of the cylinder. The inner radius of the cylinder is R = 100mm and 
the distance between the inner face of the impellers is H = 180mm, which gives an aspect ratio of H/R = 1.8. The 
turbulence properties (anisotropy, fluctuations, dissipation) are influenced by the geometry of the impellers, i.e their 
non dimensional radius Rt, the oriented angle a between the blades (see Fig. [^, the number n of blades and their 
heights hb US]. In the present paper, we consider only impellers with hb/R = 0.2 and a = ±72°. Those impellers 
are the so-called “TM60” (with 16 blades) and “TP87” impellers (with 8 blades), the characteristics of which are 
summarized in Table ^ They are essentially similar, except for their number of blades and the material they are made 
of. A single impeller can be used to propel the fluid in two opposite directions, respectively associated to the concave 
or convex face of the blades going forward. This can be taken into account by a change of sign of the parameter a. 
In the sequel, we denote (—) {resp. (±) ) an impeller used with the concave {resp. convex) face of its blades going 
forward. The impellers are driven by two independent motors which can rotate at frequencies up to lOHz. The motor 
frequencies can be either set equal to get exact counter-rotating regime, or set to different values /i 7 ^/ 2 . To change 
the viscosity, we have used either water or glycerol at different dilution rates. 


B. Control parameters 


In the sequel, we choose R and 11 ^ = ( 7 r(/i ± / 2 )) ^ as unit of length and time. The von Karman experiment is then 
characterized by two control parameters: 

• the Reynolds number: 

Re = 7r(/i ± f2)R^v~^, 

where v is the fluid kinematic viscosity, ranges from 10^ to more than 10® so that we can span a full range of 
regimes, from the purely laminar to the fully turbulent one. 

• the rotation number: 


e = 


/1-/2 

/l±/2’ 


measures the relative influence of global rotation over a typical turbulent shear frequency. Indeed, the exact 
counter-rotating regime corresponds to 0 = 0. For a nonzero rotation number, our experimental system is 
similar, within lateral boundaries, to an exact counter-rotating experiment at frequency / = (/i ± / 2 )/ 2 , with 
an overall global rotation at frequency (/i — /2)/2 |15[ [22] . In our experiments we vary 6 from —1 to ±1, 
exploring a regime of relatively weak rotation to shear ratio. 

Table ]^ summarizes the parameter space explored in this paper. The Reynolds variation is done at 0 = 0 while the 
rotation variation has been explored at Re « 10®. 
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impellers 

material 

number of blades 

a (in degrees) 

Re 

e 

TM60(+) 

stainless steel 

16 

72 

[10^10®] 

0 

TP87(+) 

polycarbonate 

8 

72 

10® 

[-0.5, 0.5] 

TP87(-) 

polycarbonate 

8 

-72 

10® 

[-0.5, 0.5] 


TABLE I: Parameter space explored in this paper. 


C. Measurements 

The set-up allows for both global and local flow diagnostics. Torque (global) measurements at each impeller are 
performed with SCAIME technology and provide values over the kHz range of Ci and C 2 , being respectively the 
torque applied to the bottom and top shafts. Following the procedure described in m, they are calibrated using 
measurements at different mean frequencies, so as to remove spurious contributions from genuine offsets or mechanical 
frictions. From this, we compute the nondimensional value Kpi and Kp 2 of the torque as: Kpi = where 

p is the density of the working fluid. 

Local measurements of the velocity held of the flow have been made by PIV techniques in the stationary regime. The 
typical size of the particles used is a few tens of micrometers and their density is 1.4. Two cameras take between 
600 and a few thousand successive pictures of the flow at a 15Hz frequency. The resolution of our camera frame is 
1600x1200 pixels, and the reconstruction is done using peak correlation performed over overlapping windows of size 
16 to 32 pixels. As a result, we get measurements of velocity held on a grid of approximate size 60^ in a vertical plane 
containing the axis of symmetry (Oz), in a cylindrical system of coordinate. The maximum spatial resolution we can 
reach for the velocity held with this set-up is therefore of the order of 200/60 « 3mm, about 10 to 100 times larger 
than the dissipative scale. For more details about the experimental set-up or the measurement techniques see for 
instance m- We can therefore estimate the derivatives of the velocity held only along the r and z directions, but we 
do not have access to derivatives along 6. As a consequence, we will either set these derivatives to zero in the sequel, 
or take them into account using hypothesis of statistical axisymmetry. We have checked that the two procedure give 
essentially the same result. In any case, we use incompressibility to estimate dgug. 


D. Diagnostics 

1. Flow geometry 

We are going to study different types of von Karman flows (see Fig. [^. These types of flow happen depending 
on the forcing condition (+) or (-), and whether the system undergoes a spontaneous phase transition (bifurcation) 
or not m- There is, then, four types of flows: one corresponding to the (+) forcing condition where a phase 
transition cannot be observed (Fig. Ha)), one corresponding to the (-) forcing condition where a phase transition is 
not observed (Fig. [^b)), and two more corresponding to the two states of the flow that can be observed once the flow 
has undergone its phase transition (Fig. Hc)-(d)). The difference between these flows can be characterized through 
our PIV measurements by their mean velocity profile in the vertical plane of measurements (see Fig. H. In the first 
two types, the mean flow is symmetric with respect to the equatorial plane z = 0 and there is a strong shear layer in 
the middle. In the bifurcated states, the flow is no longer symmetric with respect to the equatorial plane. It consists in 
a one-cell flow in the vertical direction, with a strong shear layer at the impeller that rotates in the direction opposite 
to the orthoradial mean flow. The two bifurcated states are symmetric to each other. The interest of considering 
these different mean flow geometries is that they are characterized by well resolved (resp. badly resolved) shear layer 
for the symmetric state (resp. bifurcated state). Since we expect an important fraction of energy dissipation to be 
localized where there is a strong velocity gradient, this difference may be a large source of error in local estimates 
based on PIV measurements, as we demonstrate later. 


2. Global dissipation 

Using torque measurements, we can get an accurate estimate of the global power injected into the flow. Indeed, since 
we study a stationary situation, this input must balance the rate of energy dissipation within the flow. This has 
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FIG. 2: The four types of flow characterized by their mean velocity profiles. Arrows represent the velocity field in the radial 
plane while colours represent the orthogonal component to that plane, a) The flow is forced under the (+) condition, no 
phase transition can be observed, b), c) and d) show the three different states of a flow forced under the (-) condition, b) no 
phase transition has occurred, c) (resp. d)) the shear layer has been sent downwards (resp. upwards). In each of these flows 
Re « 3 X 10®, and the two TP87 impellers rotate in opposite directions at 6 — 0. 


been checked in a scale 4:1 version of our experiment (see m)- From these measurements, we get the global mean 
dissipation rate as D = 27r(C'i/i + C 2 / 2 ). From this, we compute diagnostics: 

• the dimensionless mean dissipation rate: 

• the dimensionless mean dissipated power per unit mass: 



VR 


(15) 


which should not be confused with the Levi-Civita symbol. V — ttR^H is the volume of the experiment. 

These quantities depend on the Reynolds number, the rotation number, the characteristics of the impellers and the 
mean flow geometry |23| . Examples of variation of e as a function of Re are provided in Fig. I^a) at 0 = 0, with 


























































TP87 and TM60 impellers, for the different flow geometries illustrated in Fig. At low Reynolds numbers, the 
dissipation rate decays as Re~^, until Re « 300 where the turbulence sets in. The dissipation rate then levels off at 
a value which depends upon the flow geometry: it is the lowest for the symmetric (+) flow, then increases for the (-) 
symmetric flow and is the largest for the two bifurcated (-) flows. In Fig. [^b), the mean dissipation rate is shown as 
a function of 6 for TP87 (±) impellers at Re = 3 x 10®. In the (+) case (pink disks), it may be seen that the minimum 
of e is obtained at 9 — 0. When 9 is varied from, say, 0 to 0.5, e increases and reaches a value twice as large as its 
minimum. In the (—) case, there is a discontinuity of the energy dissipation at 9 = Q due to the global symmetry 
breaking. We see that the symmetric state (blue rhombi) dissipates less energy than the two bifurcated states (red 
triangles). Another striking feature of Fig. ib), is the coexistence of two branches of dissipation, corresponding to 
the two bifurcated states, for a certain range of values of 9. This coexistence vanishes around 9 = ±0.2 (this value 
depends on the Reynolds number, see ED)- 


ip) 



Re e 


FIG. 3: Plots of the dimensionless dissipated power per unit mass: a) As a function of Re for TP87 and TM60 impellers for 
the four geometries depicted in Fig. [^at 9 = 0. Dotted line: Fit e ~ 37Re~^. Pink disks: symmetric state (±). Blue rhombi: 
symmetric state (-). Red triangles: bifurcated states (-). b) As a function of 9 for TP87 impellers. Re = 3 x 10®. Pink disks: 
symmetric state (+). Blue rhombi: symmetric state (-). Red triangles: bifurcated states (-). 


3. Local production and dissipation rates 

Since the von Karman flow is neither isotropic, nor homogeneous, it is interesting to study more locally the production 
and dissipation of energy. Our set-up being statistically axisymmetric, it is sufficient to perform the study in the 
vertical plane including the rotation axis, where we perform our velocity measurements. We thus consider the main 
stationary terms of the local energy balance Eq. (|^, and divide them by the experimental volume V = ttR^H to get 
three local quantities: 

• the local production rate per unit mass: 


r = -^d,jR (16) 

Integrating over the whole volume, one may then get from this the corresponding total production rate per unit 
mass Ft: 


Tt 


1 

V 


diJ'-dV 



( 17 ) 


where the last equality comes from Green-Ostrogradsky formula. 






9 


• the positive local dissipation (transfer) rate per unit mass: 

n = y + 2vSijS^i^ . (18) 

The first term represents the energy dissipation due to energy transfers towards SGS while the second one is 
the laminar energy dissipation. 

• the positive local singularity dissipation (transfer) rate per unit mass: 

Ifrifi = (19) 

As before, we denote by Hr and TioRT the total dissipation rate per unit mass, obtained by volume integration 
of n and 

With our stereoscopic PIV measurements, we have access to all terms to compute these quantities, except for the 
pressure P in J*, and the terms involving derivatives with respect to the azimuthal angle. In the sequel, we will 
present maps of the local quantities without the terms we do not have access to. In other words, we neglect all terms 
involving pressure and azimuthal derivatives. Moreover, the computation of the transfer terms requires the calibration 
of the constant C involved in the gradient model (see Eq. @). In section IV A, we describe a calibration procedure 
using the angular momentum balance equation. 


IV. DATA PROCESSING 


A. Calibration of the gradient model using angular momentum balance 


The first step is then to find the value of the C constant in order to be able to estimate the contribution of the 
terms containing in (161 and (18l. For this, we follow the work of [53] who has shown, by using high-resolution 


Laser Doppler Velocimetry (LDV) measurements, that in a symmetric situation (0 = 0), the vertical flux of angular 
momentum is a constant equal to the torque injected by the impellers. This result stems from the z-component of 
the angular momentum balance Eq. ([^ integrated over a volume V(z) that describes a cylinder extending from the 
bottom impeller, to an altitude h{z) = H/2 + z (see Fig. |^. In this case, the global z-component of the angular 
momentum balance equation ([^ reads 


Kpi = —Re ^ I / {ue + rdrUg)dS — / rdzU0dTi \ + / r{urUeT'^^)dS + / r(ugUz + T^^)dY^ 
\JS(z) J Js{z) Jt.{z) 

= -Ky{z)+ [ r{ugUz + T^^)dY., (20) 

Js{z) 

where Il(z) is a surface at altitude h{z) from the bottom impeller, and S{z) is the lateral boundary (see Fig. |^. As 
discussed in |24| , this equation states that the angular momentum transmitted by the motor to the fluid (measured by 
the dimensionless torque at the bottom Kpi) is either evacuated through the lateral boundary, or transported to the 
upper layers of the fluid to be received by the upper motor as a drag. It was further shown in |24| that in a symmetric 
situation (0 = 0), where the torques at the bottom and at the top are equal, the Ky{z) term is negligible, so that 
there is a constant flux of angular momentum from the bottom to the top. For each 0 = 0 case, we thus compute the 
quantity 


[ r{uguz +T^^)dS, ( 21 ) 

and adjust the constant C, so that $i^(z) = Kpi = Kp2 = Kp for any z. In Fig. we present the results of our 
computations for TM60(+) impellers. The adjustment of the constant C is done statistically, by assuming that it 
does not depend on the Reynolds number (as it is supposed to depend only on the velocity reconstruction algorithm). 
We then fix = A^Az, where = 2i?/58 and Az = {H — 2/ib)/63 being the spatial resolutions of our PIV and we 
adjust C using 38 symmetrical flows at different Reynolds numbers so that their statistical mean provides a constant 
value of 4>(z)/A'p = 1 within the (statistical) standard deviation. 
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FIG. 4: Left: Control volume for the angular momentum budget. Right: Plot of the normalized vertical flux of vertical angular 
momentum as a function of 2 . The dark grey area represents convective transport due to the mean flow while the light grey 
area represents convective transport due to fluctuations. Filled squares represent the sum of the two contributions so that the 
total convective transport is almost constant with 2 and equals to Kp. Each point has been obtained by taking the mean of 
several computations obtained from several flows at different Re. The error bars represent the standard deviation. 


The optimal value of C has to be taken equal to C « 4. This is around 50 times more than the conventional choice 
made in [12 ng for a Gaussian filter. We have identified several factors that may explain this difference: our data are 
not filtered in a Gaussian way; since we ignored azimnthal derivatives, the constant must be higher to compensate; 
the tnrbnlence is neither isotropic nor homogenons. 

Dividing the flux into the contribntion due to the time-average flow and the contribution dne to the fluctnations, 
we also recover that near the impellers the flux is only due to the mean flow while at the center it is almost only 
dne to fluctnations [24]. We also checked that the diffnsive terms in (201, going as Re~^, are negligible at high 
Reynolds number. Likewise, the flux through S{z) is also small, meaning that almost all of the flux goes in the 
vertical direction. 


B. Implementation of the Duchon-Robert formula 


As we have seen in section HE ^i{u) provides the local dissipation for a given velocity field u, at a given length scale 
Since this formnla is a generalization of the Karman-Howarth equation to any kind of flow, it provides an estimate 
of the snbgrid energy transfer as long as the considered scale I is in the inertial range. Moreover, the considered scale 
mnst be snfhciently large with respect to the PIV smallest resolved scale, so as to gnarantee statistical convergence 
through sufficient average in the scale space (on the sphere of radius t). To check these two points, we present on 
Fig. ^e{u) in two plots, averaged in the radial (resp. vertical) direction, as a function of 2 (resp. a;) and i, for an 
experiment nsing TM60(+) impellers, at very large Reynolds nnmber {Re « 8 x 10®). We see that this quantity is 
close to zero at large scales {i > 0.4), bnt that it increases at small scales in the domain |a:| < 0.4, \z\ < 0.7, i.e. at 
the location of the median shear layer. For i between 0.1 and 0.15 (i.e approximately 4 to 5 times the smallest scale 
resolved by onr PIV set-np), there is the start of a saturation, indicating the beginning of the inertial range. While 
the extent of the inertial range is likely to vary (increase) with the Reynolds number, its largest scale is likely to be 
independent of the Reynolds number, as long as the flow is turbulent. Indeed, as discussed in [^, the geometry of 
the largest scales of the von Karman flow appear fairly independent of the Reynolds nnmber, except aronnd Re = 10® 
where they may experiment abrupt changes due to the equivalent of a phase transition m- Since in the Kolomogorov 
pictnre energy cascades from large to small scales, it is reasonable to assume that the beginning of the inertial range is 
solely determined by the large scale topology, thereby becoming independent of the Reynolds nnmber (except maybe 
aronnd Re = 10®). To check this, we report on Fig. ^a) the comparison between the total dissipated power I\.drt 
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at £ = 0.1 (blue symbols) and direct measurements of the injected power (black symbols). We observe that our 
estimates are in good agreement with direct measurements, especially for symmetric flows. In Fig. ^b), the results of 
our computations to estimate the global dissipated power nsing the LES method are displayed. These resnlts will be 
discnssed in more detail later, but we can already observe that we obtain a good agreement with direct measurements. 
We can conclnde from all our computations that the method nsing the DR formnla may be seen as an interesting 
alternative to the widespread LES-PIV method, since it relies on very few arbitrary hypotheses. We explore its 
performances with other flow geometries in the sequel. 


(b) 


< ^e(u) >z 




FIG. 5: Vertical (a) and radial (b) average of the local dissipation &e{u) for TM60(+) at Re = 8 x 10® as a function of x, z 
and £. 
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FIG. 6: Comparison between direct measurements of energy injection obtained using torque measurements (black symbols) 
and PIV estimates at various Reynolds numbers, for TM60 and TP87 impellers for the four different mean state geometries 
displayed in Fig. disks: (+) symmetric. Rhombi: (-) symmetric. Up triangles: (-) shear layer sent downwards. Down 
triangle: (-) shear layer sent upwards, a) energy dissipation Hdrt using the DR method (blue symbol), b) energy injection 
Ft (red symbols) and dissipation IIt (blue symbols) using the LES-PIV method. The dotted line represents the "laminar fit" 
e = 37Re~^. The estimates are computed based on 2 to 15 realizations of the experiment where at least 600 instantaneous 
velocity snapshots have been taken for each of them. The symbols represent the mean of our computations while the error bars 
represent the standard deviation. 
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V. RESULTS 
A. Symmetric case 0 = 0 

a. Energy production and transport To understand the flux of energy in the von Karman flow, it is interesting to 
focus first on the symmetric case {9 = 0), specifically on the case when both impellers rotate in the (+) sense with 
the same frequency, in a stationary regime. In such regime, the dissipated power equals the injected power. Through 
Green-Ostrogradski theorem, the total energy production Ft is equal to the entering flux of at the boundaries. In 
this symmetric case, most of the flux is provided by the component J^, and appears to provide a fairly good estimate 
of the injected power at large enough Reynolds numbers (Re >= 10®). Indeed, we show on Fig. §b) the results of 
our computations for Tt (empty red symbols). We see that for Re >= 10®, these estimates coincide within 20%. 
However, at lower Reynolds numbers (Re « 4 x 10®), the PIV estimates only capture 55% of the actual injected 
power. This may be due to the fact that as the Reynolds number is decreased, an increasing part of the injected 
power is either through pressure effects, azimuthal variations or viscous boundary layers that are not resolved by our 
measurements. 

Given that the contribution dominates the total energy production, it is interesting to focus on the spatial variation 
of this quantity. Fig. j^a) shows such a local map of for a symmetric flow at Re « 3 x 10^. 




FIG. 7: Typical maps of a) mean energy transfer through the system, b) contribution of the mean flow, c) contribution of 
the fluctuations 

As can be seen, at the center of the cylinder, i.e. x/R G [—0.8 0.8], big structures reflect the advection of energy 
towards the impellers, through the mean Ekman pumping by the impellers. In contrast, at the walls, smaller structures 
of opposite sign are observed, reflecting the injection of energy within the flow. These two kinds of structures are 
mirrored by the recirculation cells that we observed on Fig. The local structure of J* can be further used to get 
information about energy transport in the flow. Indeed, performing an integration of equation Q over the height- 
varying volume V(z) (see Fig. [R and setting dtE = 0, we can get an equation for the energy transport in the flow 
as 


PdS, = -Kpi{e + l)-K(,{z)+ / J^dS, 


( 22 ) 


where AT' is the contribution due to the lateral boundaries at height 2 . Ignoring this contribution, we can write 


J^dSi = -Kpi + ^e(z)j 


(23) 


where ^^(z) = J^dS. This quantity is displayed on Fig. 8 . 

We observe that near the impellers, the total flux of energy ^lled black squares) equals the energy injected by the 
impellers. The change of sign comes from the fact that the impeller at the bottom injects energy in the +z direction 
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FIG. 8: Plot of the normalized vertical flux of energy as a function of z. The dark grey area represents convective transport 
due to the mean flow while the light grey area represents convective transport due to fluctuations. Filled squares represent the 
sum of the two contributions so that the total convective transport is a decreasing function of 2 . Each point has been obtained 
by taking the mean of several computations obtained from several flows, and the error bars represent the statistical standard 
deviation. 


whereas the upper impeller injects energy in the —z direction. At the center of the experiment the flux is zero, 
meaning that on average there is not any energy going from one half of the cylinder to the other through the shear 
layer. Finally, it is interesting to see that, as before, the mean flow plays an important role in the transport of energy 
near the impellers, whereas near the shear layer energy is carried by fluctuations. A slight difference is that near the 
impellers, as the mean flow sends energy towards the center of the cylinder, fluctuations create a flux which goes in 
the opposite direction and tries to send energy back to the impellers. 

Altogether, the local map of injected power is shown on Fig. |^a). We see that near the impellers, the divergence 
term brings energy into the system as is expected. We also see that energy leaves the center of the recirculation cells 
and tends to be advected towards the center of the impeller through Ekman pumping. 

b. Energy dissipation We now turn to the total energy dissipation estimated through the LES method IIt’. 

Its values at various Reynolds numbers in the symmetric case are reported in Fig. §b). We observe that it is in 
very good agreement with direct measurements at Re « 4 x 10^. At such Reynolds number, the dissipative scale is of 
the order of 1mm, close to the PIV resolution. This result is therefore in agreement with the observation of Tokgoz 
et al. |T] obtained for a Taylor-Couette flow. At larger Reynolds numbers, the estimates using the LES method are 
in good agreement with respect to direct measurements since we are able to capture up to 90% of the actual energy 
dissipation. The map of the local dissipation with this method is provided in Fig. ib), and appears fairly uniform 
across the vessel, with larger intensity along the vertical axis around r = 0 and near the lateral boundaries around 
z = ±0.3. _ 

This map can be compared with the one obtained using the DR formula at ^ = 0.1 Fig. [^d). One observes 

the same localized structures of energy dissipation near the walls at x = ±0.6, which are symmetrically distributed 
with respect to the mid-plane. These structures are statistically significant, since they are not observed on plots of 
the instantaneous local dissipation (see Fig. |^c)). They may, therefore, trace the intense vortices of the shear 

layer. In addition, one observes a clear localization of energy injection (red areas) at the tip of the impellers, with 
an energy dissipation in the middle part of the cells. Regarding instantaneous maps, it is interesting to note that it 
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FIG. 9: Top: Maps of the energy production rate (a) and energy transfer rate (b) for the same flow as in Fig. [^a), using 
LES-PIV estimates. Bottom: Instantaneous (c) and time averaged (d) maps of the local dissipation &t{u) at ^ = 0.1 for the 
same flow. Areas where energy accumulates are represented in red, while those where energy leaves are represented in blue. 


is characterized by intense, localized events, that can reach 10 to 20 times the maximum average local dissipation. 
Whether these intense events are connected to localized quasi singular structures is an interesting open question that 
we leave for future work. 

Overall, these dissipation maps provide clear evidence that the maximum energy dissipation lies within the shear layer 
and that the DR formula provides a better estimate of the energy dissipation than the LES method in the symmetric 
case. 

c. Description of the energy cycle Altogether, our results regarding energy production, transport and dissipation 
can be summarized into a simple picture of the "energy cycle" for the von Karman flow, which is sketched in Fig. 
10 the energy is advected to the impellers via the Ekman pumping, the flow is then accelerated inside the impellers 
and expelled at the walls, providing an energy injection towards the mixing layer. In that mixing layer, turbulent 



























15 


fluctuations dissipate an energy equal to e. The fluid is then pumped again into the impellers for further reinjection, 
closing the energy cycle and providing a stationary situation with energy injection and dissipation equal to e. 

In the sequel, we study the influence the flow topology and the forcing conditions onto this energy cycle. Since in 
this picture most of the energy is dissipated within the middle shear layer, there is good hope that we can capture 
its main contribution by the PIV measurements, provided that the shear layer is not too close to the impeller since 
in that case, the PIV measurements cannot resolve the flow. We check in section V B-C that this is indeed the case, 
and provide detailed informations about the localization of energy dissipation to complete the picture of the energy 
cycle. 



FIG. 10: Energy cycle in the von Karman experiment: the energy is advected to the impellers via the Ekman pumping, the 
flow is then accelerated inside the impellers and expelled at the walls, providing an energy injection e towards the mixing layer. 
In that mixing layer, turbulent fluctuations dissipate an energy equal to e. The fluid is then pumped again into the impellers 
for further reinjection, closing the energy cycle and providing a stationary situation with energy dissipation e. 


B. Influence of the flow topology 


Let us now consider the case when both impellers rotate in the (-) sense at 0 = 0. For this type of forcing, at sufficiently 
high Reynolds numbers, there is coexistence of three different flow geometries for TP87 type impellers. 

Considering the case where the flow has not undergone any phase transition. Fig. [^shows that for this type of forcing 
we measure a dissipation almost three times bigger than what was observed with the previous forcing condition. In 
this symmetric case, our estimates of the injected and dissipated power are within 20% of the measured value using 
the LES-PIV method, whereas we reach 98% of the actual dissipation rate of energy with the DR formula (see Fig. 
[^a)). Local maps of injected and dissipated power are plotted on Fig. They correspond to the flows (b) displayed 
on Fig. 1^ where the forcing is in the (-) sense. Here Fig. (a) and (b) represent LES-PIV estimates while Fig. 11 
(c) and (d) represents instantaneous and time averaged maps using the DR method. Here, the remarks are essentially 
the same as in Fig. near the impellers, the divergence term brings energy into the system while energy leaves the 
center of the recirculation cells to accumulate at the center of the cylinder. The dissipation term does not change 
much either and is approximately constant throughout the plane of measurements. A noticeable difference however, 
is that the color scale is bigger for this configuration. We thus recover the fact that the flow dissipates more energy 
than for a forcing in the (+) sense and, as a consequence, the torque of the impellers has to be higher. 

The situation changes for the bifurcated (-) states. The shear layer is now very close to the upper or lower impeller, 
depending on the state "chosen" by the system. The mean velocity held for these two states is essentially the same, 
differing only by the transformation z —> —z. The measured injected power is also the same in these two states. 
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FIG. 11: Top: Maps of the energy production rate (a) and energy transfer rate (b) for symmetric (-) geometry corresponding to 
the flow in Fig. Sb), using LES-PIV estimates. Bottom: Instantaneous (c) and time averaged (d) maps of the local dissipation 
&i{u) at ^ = 0.1 for the same flow. Areas where energy accumulates are represented in red, while those where energy leaves 
are represented in blue. 


However, in terms of PIV-estimated injected or dissipated power, a clear asymmetry occurs. Indeed, when the shear 
layer is sent downwards (up-pointing triangles), the LES-PIV method provides good estimates of both the injected 
and the dissipated power (see Fig. ®b)), while the DR-method provides about 65% of the dissipated power. In the 
case where the shear layer is sent upwards (down pointing triangles), however, both the LES-PIV method and the 
DR method totally fail to reproduce the measured dissipated energy (giving the wrong sign). Since there is not any 
asymmetry observed in the mean flow, this difference must be attributed only to fluctuations and tiny asymmetries 
of the experimental set-up (laser sheet location, localization of cameras, ..). This points out the the importance of 
resolving the shear layer in the PIV-estimates of the injection or dissipation energy rate. We thus conclude that 
neither the LES-PIV, nor the DR method are appropriate in cases when the shear layer is at the location of the 
impellers, where it cannot be resolved by our PIV set-up. 

Local maps of injected and dissipated power in the bifurcated states are plotted on Fig. They correspond to the 
flows (c) and (d) displayed on Fig. where the forcing is in the (-) sense. Maps (a) (c) and (e) correspond to a flow 
where the shear layer has been sent near the lower impeller. As a consequence, we observe that energy enters from 
the top and is advected downwards to the shear layer. The map of the dissipation term for this flow looks quite like 
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FIG. 12: Maps of the energy production rate and energy transfer rate for the two bifurcated geometries (-): Left: for the flow 
represented in Fig. ^c). Right: for the flow represented in Fig. [^d). Top line represents the estimated injected energy, using 
LES-PIV method. Middle line represents the estimated dissipated energy using the LES-PIV method. Bottom line represents 
the estimated dissipated energy using the DR method. Areas where energy accumulates are represented in red, and those where 
energy leaves are represented in blue. 
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the production term. We see that most of the dissipation we are able to capture happens at the center of the cylinder 
but, as we said before, we miss all the dissipation that happens near the lower impeller where there is the shear layer. 
At first sight, maps (b) (d) and (f) are just symmetric of maps (a) (c) and (e) with respect to z —>■ —z. However, by 
looking closely, tiny differences are observed, especially near the walls and impellers. These differences explain the 
asymmetry between the two bifurcated states shown on Fig. 


C. Influence of the forcing asymmetry 6 

We now investigate the case where the impellers are not exactly counter-rotating, so that the parameter 0 varies from 
—0.5 to 0.5. In the (-) case, the geometry observed at |0| > 0 is always bifurcated (see Fig. [^b)), with the shear layer 
located at the top or bottom impeller. In the (+) case, the transition is more gradual: the shear layer is increasingly 
shifted upwards (resp. downwards) as 9 goes from 0 to 1 (resp. —1), allowing for finer tests of the accuracy of our 
PIV-estimates as a function of the flow geometry. 
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FIG. 13: Comparison between direct measurements of energy injection obtained through torque measurements (black symbols) 
and PIV based estimates at various 6, for TP87 impellers for the four different mean state geometries of Fig. 2: disks (+). 
Rhombi: (-) symmetric. Up triangles: (-) shear layer sent downwards. Down triangle: (-) shear layer sent upwards, a) energy 
dissipation using DR method (blue symbols), b) energy injection Ft (red symbols) and dissipation IIt (blue symbols) using 
the LES-PIV method. The estimates are computed based on 1 to 2 realizations of the same experiment where at least 600 
instantaneous velocity snapshots have been taken for each of them. The symbols represent the mean of our computations while 
the error bars represent the standard deviation. 


On Fig. 13 are plotted our estimates for the total injected (red symbols) and dissipated power (blue symbols) using 
both the LES-PIV method (Fig. [T^b)) and the DR formula (Fig. [T^a)). For each 9, our computations are based on 
sets of at least 600 instantaneous velocity snapshots. Because of the symmetry z —>■ — z, we expect all estimates to be 
symmetric with respect to 0 —>■ —9, provided the statistics are well converged and that the shear layer is sufficiently 
resolved. We see that this is indeed the case for (+) forcing condition, where the shear layer always lies in between 
the two impellers. However, we observe that both methods become inaccurate when |0| is too high i.e \9\ > 0.3. 
For the (-) forcing condition, the estimates give very good agreement for the symmetric state. However, when the 
phase transition occurs, one observes the same asymmetry as in the 0 = 0 case between the two possible states: when 
the shear layer is sent downwards we get a good agreement between measurements and PIV estimates with the LES 
method. However, the DR method systematically underestimates the dissipation. In the case where with the shear 
layer is sent upwards, the estimates are really bad with both methods. 

Corresponding maps of the injected and dissipated power are provided in Fig. 
sense. They obey qualitatively the same behavior as in the symmetric case, so that the energy cycle description is 
qualitatively the same, with the shear layer location being moved as 0 varies. 


14 at 0 = —0.1 in the (+) and (-) 
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FIG. 14: Maps of the energy production rate and energy transfer rate at 9 = —0.1 for the two rotation directions: Left: (+). 
Right: (-). Top line represents the estimated injected energy, using LES-PIV method. Middle line represents the estimated 
dissipated energy using the LES-PIV method. Bottom line represents the estimated dissipated energy using the DR method. 
Areas where energy accumulates are represented in red, and those where energy leaves are represented in blue. 
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VI. CONCLUSION 

In this paper, we investigate the energy cycle of a turbulent von Karman flow. In this kind of set-up, it is possible 
to control the mean energy flux inside the flow in order to get a statistically stationary regime. It is also possible to 
monitor the torque and frequency of the impellers. In our case, we use Particle Image Velocimetry (PIV) measurements 
to study local and global energy transfers inside the flow. This kind of measurements gives us access to the effective 
velocity field on a grid that does not resolve the dissipative scale. Therefore, we use a Large Eddy Simulation (LES) 
approximation to model the influence of unresolved scales. This procedure involves a free parameter which has to 
be calibrated for our set-up. This is achieved by imposing angular momentum balance at the smallest resolved 
scale [23]. After deriving an energy balance equation at a fixed scale £, we proceed to estimate four quantities from 
our PIV measurements: the local and global mean power injected by the impellers and the local and global mean 
dissipated power. This computation is performed for various Reynolds numbers and for various flow topologies. 
These PIV-estimates are then compared with direct injected power estimates provided by torque measurements 
at the impellers. The agreement between PIV-estimates and direct measurements depends on the flow topology. In 
symmetric situations, we capture up to 90% of the actual energy dissipation. However, our results become increasingly 
inaccurate as the shear layer responsible for most of the dissipation approaches one of the impellers, and cannot be 
resolved by our PIV set-up. At the same time, we explore a new method for PIV-estimates of the energy dissipation, 
based on the work of Duchon and Robert [3] that generalizes the Karman-Howarth equation to nonisotropic, non 
homogeneous flows. This method provides parameter-free estimates of the energy dissipation as long as the smallest 
resolved scale lies in the inertial range and the shear layer is resolved by the PIV set-up. With this method, we obtain 
a very good agreement between PIV estimates and direct measurements, and we are able to capture up to 98% of 
the actual dissipated power in symmetric situations. However, this method also gives increasingly inaccurate results 
as the mixing layer approaches one of the impellers. These results are used to evidence a well-defined stationary 
energy cycle within the flow in which the energy is injected by the top and bottom impellers towards the shear layer. 
There, turbulent fluctuations dissipate energy, and the flow is then pumped towards the impellers, closing the energy 
cycle. 
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